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We study heavy hadron spectroscopy near heavy meson thresholds. 
We employ heavy pseudoscalar meson P and heavy vector meson P* 
as effective degrees of freedom and consider meson exchange potentials 
between them. All possible composite states which can be constructed 
from the P and P* mesons are studied up to the total angular momen- 
tum J < 2. We consider, as exotic states, isosinglet states with exotic 
J PC quantum numbers and isotriplet states. We solve numerically the 
Schrodinger equation with channel-couplings for each state. We found 
B^B^ molecule states for I G (J PC ) = l + (l + ~) correspond to the masses 
of twin resonances Z b (10610) and Z b (10650). We predict several possible 
;g(*)p/ -* bound and/or resonant states in other channels. On the other 
hand, there are no D ( *- ) D ( '*' ) bound and/or resonant states whose quantum 
numbers are exotic. 
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1 Introduction 



Recently, the exotic twin resonances Zf,(10610) and Zb(10650) are discovered in the 
processes T(55) -> 7T7rT(n£)(n = 1,2,3) and T(55) ->■ Trnh b (kP)(k = 1,2) by 
Belle group [TJ [2]. The reported masses and widths of the two resonances are 
M(Z b (10610)) = 10607.2±2.0MeV, r(Z b (10610)) = 18.4±2.4 MeV and M(Z b (10650)) = 
10652.2±1.5 MeV, r(Z b (10650)) = (11.5±2.2) MeV. These resonances have some in- 
teresting properties. First of all, ZbS have exotic quantum numbers I G (J P ) = 1 + (1 + ). 
Since Z b 's are isotriplet states, they need four quarks as minimal constituents. So 
Z b 's are "genuinely" exotic states. Secondly, Z b s have exotic decay ratios. In gen- 
eral, a decay process with hf, should be suppressed in heavy quark mass limit, be- 
cause these processes require a heavy quark spin flip. Nevertheless, the decay rates 
of T(5S f ) — > Z\jix — > inrT(nS) are comparable to those of T(5S) — > Z b n — > nirhb- 
Thirdly, Z^s are "exotic twin" resonances. Their mass splitting is only 45 MeV 
(Amz b ~ 45 MeV). This scale is not typically found in usual quarkonium. These facts 
strongly suggest that Z^s have a molecular type structures as noted in Ref. [3j H] 

We can naively expect the existence of molecular states in heavy quark sectors. 
We have two reasons. One is that the kinetic term of Hamiltonian is suppressed. 
Because the reduced mass is larger in heavy mesons. Second is that P and P* are 
degenerate thanks to heavy quark symmetry because the interaction of heavy quark 
spin is suppressed in heavy quark sector. This leads the effects of channel-couplings 
become larger. 

In this paper we study pMp^ molecular states in terms of the potential model. 
We completely take into account the degeneracy of P and P* due to the heavy quark 
symmetry, and fully consider channel couplings of PW and P^* . And we consider not 
only bound states but also resonant states. 

This paper is organized as follows. In Sec. 2, we introduce (i) the tt exchange 

potential and (ii) the npcu potential between P(*) and P^* mesons. To obtain the 
potentials, we respect the heavy quark symmetry for the PWp0O7T, p(*)p(*)p an d 
pWpWo; vertices. We classify all the possible states composed by a pair of PW and 

P * mesons with exotic quantum numbers I G (J PC ). In Sec. 3, we solve numeri- 
cally the Schrodinger equations with channel-couplings and discuss the bound and/or 

t v ,,w .„«„.,...,„.. . w. , u , ^*>B and -D'^-D systems. In Sec. 4 is devoted to summary. 
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2 Method 



We employ the effective Lagrangians based on heavy quark and chiral symmetries. 
They give the interaction Lagrangians of ttPP* and 7rP*P* with 

ArPP* = 2f(PlP b * fl + P:lP b )d»7t ab , (1) 
Jn 

Cp*p* = 2i^e^v a P*JP b \d u n ah , (2) 

where /„- = 135 MeV is the pion decay constant. The coupling constant \g\ = 0.59 
is determined with reference to the observed decay width T = 96 keV for D* — > D71", 
assuming that the charm quark is sufficiently heavy. And the interaction Lagrangians 
of vPP, vPP* and vP*P* (v = p, co) are given by 

CvPP = -V2f3g v P b P^v ■ p ba , (3) 
C vPP * = -2V2\g v v,e^(plP;p-P:lP b )d v (p a ) ba , (4) 

£ vP *p» = V2(3g v P b *P^vp ba 

+i2V2\g v P:lP; v {d^n ba - 0"(p") fa ) , (5) 

where gy = 5.8 is the coupling constant for p — > tttt decay. The coupling constants 
are fixed as /3 = 0.9 and A = 0.56 GeV, which are determined by the radiative decays 
of D* meson and semileptonic decays of B meson with vector meson dominance by 
following Ref. [5]. Due to the G-parity, the signs of vertices for vPP, vPP and 
vP P are opposite to those of vPP, vPP* and vP*P*, respectively, for v — 00, while 
they are the same for v = p. 

Nex, we classify all the possible quantum numbers I G (J PC ) with isospin /, G- 
parity, total angular momentum J, parity P and charge conjugation C for the states 

which can be composed by a pair of P^*^ and mesons. The charge conjugation C 
is defined for / = or I z = components for 1 = 1, and is related to the G-parity by 
G = (— 1) 7 C. In the present discussion, we restrict upper limit of the total angular 
momentum as J < 2, because too higher angular momentum will be disfavored to 
form bound or resonant states. To derive a potential between P and P*, The pMp * 
components in the wave functions for various 

jPC 

are listed in Table [TJ We use the 
notation 2S+1 Lj to denote the total spin S and relative angular momentum L of the 

two body states of P^*^ and P^ mesons. The J PC = + ~ state cannot be generated 
by a combination of P(*) and mesons. For 1 = 0, there are many B^B*'*' 1 states 
whose quantum number J PC are the same as those of the quarkonia as shown in 
the third row of / = 0. In the present study, however, we do not consider these 
states, because we have not yet included mixing terms between the quarkonia and 
the p(*)p *^ states. This problem will be left as future works. Therefore, for / = 0, we 
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consider only the exotic quantum numbers J PC = , 1 h and 2 + ~. The states of 
1 = 1 are clearly not accessible by quarkonia. We investigate all possible J PC states 
listed in Table CD 

We obtain the potentials with channel-couplings for each quantum number I G (J PC ) 
in terms of the interaction Lagrangians. For each state, the Hamiltonian is given as 
a sum of the kinetic energy and the potential with channel-couplings in a form of a 
matrix as 



Breaking of the heavy quark symmetry is taken into account by mass difference be- 
tween P and P* mesons in the kinetic term. The explicit forms of the Hamiltonian for 
each I G (J PG ) are presented in Ref. [Bj. For example, the J PC = l + ~ state has four 



and hence it gives a potential in the form of 4 x 4 matrix. 

3 Numerical results 

To obtain the solutions of the pWp * states, we solve numerically the Schrodinger 
equations which are second-order differential equations with channel-couplings. As 
numerics, the renormalized Numerov method developed in Ref. [7] is adopted. The 
resonant states are found from the phase shift 5 as a function of the scattering energy 
E. The resonance position E r is defined by an inflection point of the phase shift 
6(E) and the resonance width by r r = 2/ (d5 / dE)E=E r - To check consistency of our 
method with others, we also use the complex scaling method (CSM) [§]. We obtain 
an agreement in results between the renormalized Nemerov method and the CSM. 

In Table [21 we summarize the result of the obtained bound and resonant states, and 
their possible decay modes to quarkonium and light flavor meson. For decay modes, 
the p meson can be either real or virtual depending on the mass of the decaying 
particle, depending on the resonance energy which is either sufficient or not to emit 
the real state of p or u meson. p*(oo*) indicates that it is a virtual state in radiative 
decays assuming the vector meson dominance. We show the mass spectrum of these 
states in Fig [TJ 

Let us see the states of isospin 1 = 1. Interestingly, having the present potential 
we find the twin states in the I G (J PC ) = l + (l + ~) near the BB and B*B thresholds; 
a bound state slightly below the BB threshold, and a resonant state slightly above 
the B*B threshold. The binding energy is 8.5 MeV, and the resonance energy and 
decay width are 50.4 MeV and 15.1 MeV, respectively, from the BB threshold. The 
twin states are obtained when the npuj potential is used. We interpret them as 
the Zb(10610) and Zb(10650) observed in the Belle experiment [H[2]. It should be 




(6) 




PP*-P*P) ( 3 £>i), P*P*( 3 Si), P*P*('lDi) 
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Table 1: Various components of the p(*)p * states for several J PC (J < 2). The 
exotic quantum numbers which cannot be assigned to charmonia or bottomonia QQ 
are indicated by y/. The + ~ state cannot be neither bottomonium nor p(*)p^ states. 
Examples of bottomonia are shown in non-exotic quantum numbers. 





components 


exoticness 






1 = 


I = 1 


o+- 




V 


V 


0++ 


PP(%), P*P*(%), P*P*{ 5 D ) 


Xbo 


V 


0~ 


^ (PP* + P*P) ( 3 P ) 


V 


V 


o-+ 


_^ (pp* _ p*p^ ( 3p o)j p*p*( 3 P ) 


Vb 


V 


1+- 


^ (pF - P*P~) ( 3 5x), ^ (PP* - P*P) (3D!), P*P*( 3 Si), P*P* ( 3 D 1 ) 


hb 


V 


1++ 


1_ (pp* + p*p) ( 3 Sl)) _^ (pp* + p*p) (3£)l)) P *p*( 5 D!) 


Xbl 


V 


1~ 


PP( 1 P 1 ), -1= (PP* +P*P~) ( 3 Pl), P-P*^), P*P*( 5 Pi), P*P*( 5 Pi) 


T 


V 


1-+ 


^ (PP* - P*P~) ( 3 Px), P*P*( 3 P!) 


V 


V 


2+- 


_^ (pp* _ P *p) (3£) 2 ) j P *p*( 3 P 2 ) 


V 


V 


2++ 


PP( 1 P 2 ), ^3 (PP* +P*P) ( 3 D 2 ), P*P*( 1 D 2 ), P*P*( 5 S 2 ), P*P*( 5 D 2 ), P*P*( 5 G 2 ) 


Xb2 


V 


2-+ 


X_ (pp* - P*P~) ( 3 P 2 ), ^ (pF - P*P~) ( 3 P 2 ), P*F( 3 P 2 ), P*P*( 3 P 2 ) 


Vb2 


V 


2~ 


^ (pp* + p*p) ( 3 P 2 ), ^ (pp* + p*p) ( 3 P 2 ), P*P*( 5 P 2 ), P*P*( 5 P 2 ) 


1pb2 


V 
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emphasized that the interaction in the present study has been determined in the 
previous works without knowing the experimental data of Z^'s. 

Several comments are in order. First, the bound state of lower energy has been 
obtained in the coupled channel method of BB and B*B channels. In reality, how- 
ever, they also couple to other lower channels such as tt/ij,, 7rT and so on as shown in 
Table [TJ Once these decay channels are included, the bound state will be a resonant 
state with a finite width. A qualitative discussion will be given in Section 5. Second, 
when the ti exchange potential is used, only the lower bound state is obtained but the 
resonant state is not. However, we have verified that a small change in the ir exchange 
potential generates, as well as the bound state, the corresponding resonant state also. 
Therefore, the pion dominance is working for the BB and B*B systems. Third, it 
would provide a direct evidence of these states to be BB and B*B molecules if the 
BB and B*B decays are observed in experiments. Whether the energies are below 
or above the thresholds is also checked by the observation of these decays. 

In other channels, we further predict the B^B^ bound and resonant states. The 
I (J PC ) = 1~(0 ++ ) state is a bound state with binding energy 6.5 MeV from the 
BB threshold for the ir exchange potential, while no structure for the npu potential. 
The existence of this state depends on the details of the potential, while the states 
in the other quantum numbers are rather robust. Let us see the results for the latter 
states from the wpu potentials. For l + (0 ) and 1~(1 ++ ), we find bound states with 
binding energy 9.8 MeV and 1.9 MeV from the BB threshold, respectively. These 
bound states appear also for the tt exchange potential, though the binding energy of 
the 1 _ (1 ++ ) state becomes larger. The 1 _ (2 ++ ) state is a resonant state with the 
resonance energy 62.7 MeV and the decay width 8.4 MeV. The 1 + (1 ) states are 
twin resonances with the resonance energy 7.1 MeV and the decay width 37.4 MeV 
for the first resonance, and the resonance energy 58.6 MeV and the decay width 27.7 
MeV for the second. The resonance energies are measured from the BB threshold. 
The l + (2 ) states also form twin resonances with the resonance energy 2.0 MeV 
and the decay width 3.9 MeV for the first resonance and the resonance energy 44.1 
MeV and the decay width 2.8 MeV for the second, where the resonance energies have 
are measured from the BB threshold. 

Next we discuss the result for the states of isospin 1 = 0. In general, the interaction 
in these states are either repulsive or only weakly attractive as compared to the cases 
of I = 1. The fact that there are less channel-couplings explains less attraction partly. 
Because of this, we find only one resonant state with I G (J PC ) = + (l h ), as shown 
in Fig Q] and in Table [31 The + (l~ + ) state is a resonant state with the resonance 
energy 17.8 MeV and the decay width 30.1 MeV for the npu potential. 

In the present study, all the states appear in the threshold regions and therefore 
are all weakly bound or resonant states. The present results are consequences of 
unique features of the bottom quark sector; the large reduced mass of the B«B W 
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systems and the strong tensor force induced by the mixing of B and B* with small 
mass splitting. In fact, in the charm sector, our model does not predict any bound 
or resonant states in the region where we research numerically. Because the reduced 
mass is smaller and the mass splitting between D and D* is larger. 

Table 2: Various properties of the B^B^ bound and resonant states with possible 
I (J PC ) in I = 1. The energies E can be either pure real for bound states or complex 
for resonances. The real parts are measured from the thresholds as indicated in the second 
column. The imaginary parts are half of the decay widths of the resonances, T/2. In the 
last two columns, decay channels of a quarkonium and a light flavor meson are indicated. 
Asterisk of p* indicates that the decay occures only with a virtual p while subsequently 
transit to a real photon via vector meson dominance. 



I°(J PC ) 


threshold 


E [MeV] 


decay channels 


7r-potential 


up w-potential 


s-wave 


p-wave 


i+(o+-) 










hb + 71", Xb0,l,2 + P 


l-(0++) 


BB 


-6.5 


no 


r/b + TT, T + /9 


hb + P*, Xbl + TT 


i+(o— ) 


BET 


-9.9 


-9.8 


Xbi+P* 


Vb + P, T + vr 


i-(o-+) 


BET 


no 


no 


hb + P, XbO + TT 


T+p 


1+(1+-) 


BB* 


-7.7 


-8.5 
50.4-il5.1/2 


T + 7T 


hb+7T, Xbl+P* 


1-(1++) 


BET 


-16.7 


-1.9 


T+p 


hb+P*, XbO,l+7T 


1+(1— ) 


BB 


7.0 - i37.9/2 
58.8 - i30.0/2 


7.1 - i37.4/2 
58.6 - i27.7/2 


hb+vr, Xb0,i,2+P* 


7]b+P, T + 7T 


1-(1-+) 


BET 


no 


no 


hb + P, Xbi + vr 


r/b+7r, T+p 


l+(2+") 


BET 


no 


no 




hb + 7T, Xb0,l,2+P 


l-(2++) 


BB 


63.5 - z8.3/2 


62.7-^8.4/2 


T+p 


hb+P*, Xbl,2 + 7T 


l-(2"+) 


BET 


no 


no 


h b + p 


T+p 


l+(2~) 


BB* 


2.0-^4.1/2 
44.2 - i2.5/2 


2.0-^3.9/2 
44.1 - i2.8/2 


Xhi+P* 


Vb+P, T + 7T 



4 Summary 

In this paper, we have systematically studied the possibility of the the pMp**^ bound 
and resonant states having exotic quantum numbers I G (J PC ). These states are con- 
sisted of at least four quarks, because their quantum numbers cannot be assigned by 
the quarkonium picture and hence they are genuinely exotic states. We have con- 
structed the potential of the P ( *)p < *' ) states using the effective Lagrangian respecting 
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Table 3: The BWB W bound and resonant states with exotic I G {J PC ) in I = 0. (Same 
convention as Table II.) 



I G {J PC ) 


threshold 


E [MeV] 


decay channels 


7r-potential 


up w-potential 


s-wave 


p-wave 


o-(o-) 


BB* 


no 


no 




7] h +ui, T+rj 


0+(l"+) 


BB* 


28.6 - i91.6/2 


17.8 - i30.1/2 


hb+w*, Xhi+r] 


rjb+r], T+w 


0-(2+") 


BB* 


no 


no 




K+V, Xt>o,i,2+w 



10655 

gH=g* S 

(10650) "Z b ( 10650) 



10649 



10617 



10621 



10622 



Z b (10610) 

BB* ?____„J06Q2_ 

(10604) 10594 

10596 



10606 



10566 



BB 
(10559) 

jGgPC) i+(o-) i+(i+-) l -(!++) i+(i--) i-(2++) l+(2-) + (l" + ) 



Figure 1: The B^^B^*^ bound and resonant states with exotic I G (J PC ). The 
dots with error bars denote the position of the experimentaly observed Zb's where 
M(Z b (10610)) = 10607.2 MeV and M(Z b (10650)) = 10652.2 MeV. Solid lines are 
for our predictions for the energies of the bound and resonant states when the npu 
potential is employed. Mass values are shown in units of MeV. 
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the heavy quark symmetry. The channel mixing due to the mass degeneracy of P and 
P* plays an important role to form the p(*)p bound and/or resonant states. We 
have numerically solved the Schrodinger equation with the channel-couplings for the 
P«P W states with I G (J PC ) for J < 2. 

As a result, in I — 1, we have found that the I G (J PC ) = l + (l + ~) states have a 
bound state with binding energy 8.5 MeV, and a resonant state with the resonance 
energy 50.4 MeV and the decay width 15.1 MeV. We have successfully reproduced 
the positions of Zb(10610) and Zb(10650) observed by Belle. Therefore, the twin 

resonances of Z b 's can be interpreted as the B^B 1 '*'' molecular type states. It should 
be noted that the BB -B*B, BB -B*B and B*B-B*B mixing effects are important, 
because many structures disappear without the mixing effects. We have obtained 

the other possible B(*)B ( *' ) states in I = 1. We have found one bound state in each 
l + (0 ) and 1~(1 ++ ), one resonant state in 1~(2 ++ ) and twin resonant states in 
each 1 + (1 ) and l + (2 ). It is remarkable that another two twin resonances can 

exist in addition to the Z b 's. We have also studied the B^B * states in / = and 
found one resonant state in + (l h ). We have checked the differences between the 
results from the ir exchange potential and those from the irpu potential, and found 
that the difference is small. Therefore, the one pion exchange potential dominates 
as the interaction in the B W B bound and resonant states. We also study D^D 
molecular states, but there are no bound/resonant states. Because the reduced mass 
is smaller and the mass splitting between D and D* is larger compared with bottom 
sector. 

The T(5S) decay is a useful tool to search the B^B * states. T(5S f ) can decay 

to the B^B** "* states with l + (0 ), 1 + (1 ) and l + (2 ) by a single pion emission 
in p-wave and the state with + (l h ) by oj emission in p-wave. T(55) can also 

decay to the B^B^ states with 1"(0 ++ ), 1~(1 ++ ) and 1~(2 ++ ) by radiative decays. 
In the future, various exotic states would be observed around the thresholds from 
T(5S) decays in accelerator facilities such as Belle and also would be searched in the 
relativistic heavy ion collisions in RHIC and LHC [HI HO]- If these states are fit in our 
predictions, they will be good candidates of the B^B * molecular states. 
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